arXiv:0803.3292 



How to Distinguish Dark Energy and Modified Gravity? 

Hao WeiB 

Department of Physics and Tsinghua Center for Astrophysics, 
Tsinghua University, Beijing 100084, China 

Shuang Nan Zhang 
Department of Physics and Tsinghua Center for Astrophysics, 
Tsinghua University, Beijing 100084, China 
QQ . Key Laboratory of Particle Astrophysics, Institute of High Energy Physics, 

f^*) ■ Chinese Academy of Sciences, Beijing 100049, China 

Physics Department, University of Alabama in Huntsville, Huntsville, AL 35899, USA 

(N ' 



ABSTRACT 

The current accelerated expansion of our universe could be due to an unknown energy component 
(dark energy) or a modification to general relativity (modified gravity). In the literature, it has 

■ been proposed that combining the probes of the cosmic expansion history and growth history can 
Oh' distinguish between dark energy and modified gravity. In this work, without invoking non-trivial 

dark energy clustering, we show that the possible interaction between dark energy and dark matter 
could make the interacting dark energy model and the modified gravity model indistinguishable. An 

■ explicit example is also given. Therefore, it is required to seek some complementary probes beyond 



the ones of cosmic expansion history and growth history. 
PACS numbers: 95.36.+X, 04.50.-h, 98.80.-k 
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I. INTRODUCTION 

The current accelerated expansion of our universe p], 0, H, 0, IE H S S 0, has been one of the 
most active fields in modern cosmology. There are very strong model-independent evidences [ll[ (see also 
e-g- 03) f° r the accelerated expansion. Many cosmological models have been proposed to interpret this 
mysterious phenomenon, see e.g. pj for a recent review. 

In the flood of various cosmological models, one of the most important tasks is to distinguish between 
them. The current accelerated expansion of our universe could be due to an unknown energy component 
(dark energy) or a modification to general relativity (modified gravity) [2, [l3| . Recently, some efforts have 
focused on differentiating dark energy and modified gravity with the growth function 5(z) = 5p m / p m of 
the linear matter density contrast as a function of redshift z. Up until now, most cosmological observations 
merely probe the expansion history of our universe [H, 0, 0, 0, 0, 0, 0, 0, 0, ■ As is well known, it is 
very easy to build models which share the same cosmic expansion history by means of reconstruction 
between models. Therefore, to distinguish various models, some independent and complementary probes 
are required. Recently, it was ar gue d that the measurement of growth function 5{z) might be competent, 
see e.g. [13, Q EE ES E3, EM [13, HH 113, [13, HI, HE EE EE El, [13]- If * h e dark energy model 
and modified gravity model share the same cosmic expansion history, they might have different growth 
histories. Thus, they might be distinguished from each other. 

However, the approach mentioned above has been challenged by some authors. For instance, in 26], 
Kunz and Sapone explicitly demonstrated that the dark energy models with non-vanishing anisotropic 
stress cannot be distinguished from modified gravity models (e.g. DGP model) using growth function. 
In (2?| , Bertschinger and Zukin found that if dark energy is generalized to include both entropy and shear 
stress perturbations, and the dynamics of dark energy is unknown a priori, then modified gravity [e.g. f(R) 
theories] cannot in general be distinguished from dark energy using cosmological linear perturbations. 

Here, we investigate this issue in another way. Instead of invoking non-trivial dark energy clustering 
(e.g. non- vanishing anisotropic stress), it is of interest to see whether the interaction between dark energy 
and cold dark matter can make dark energy and modified gravity indistinguishable. In fact, since the 
nature of both dark energy and dark matter are still unknown, there is no physical argument to exclude 
the possible interaction between them. On the contrary, some observational evidences of this interaction 
have been found recently. For example, in a series of papers by Bertolami et al. [281 ]. they show that 
the Abell Cluster A586 exhibits evidence of the interaction between dark energy and dark matter, and 
they argue that this interaction might imply a violation of the equivalence principle. On the other hand, 
in [29| , Abdalla et al. found the signature of interaction between dark e nerg y and dark matter by using 
optical, X-ray and weak lensing data from 33 relaxed galaxy clusters. In [421 ]. Ichiki and Keum discussed 
the cosmological signatures of interaction between dark energy and massive neutrino (which is also a 
candidate for dark matter) by using CMB power spectra and matter power spectrum. Therefore, it is 
reasonable to consider the interaction between dark energy and dark matter in cosmology. Since dark 
energy can decay into cold dark matter (and vice versa) through interaction, both expansion history and 
growth history can be affected simultaneously, similar to the case of modified gravity. Thus, it is natural 
to ask whether the combined probes of expansion history and growth history can distinguish between 
interacting dark energy models and modified gravity models. The answer might be no. So, it is required 
to seek some complementary probes beyond the ones of cosmic expansion history and growth history. 

In this note, we propose a general approach in Sec. [Til to build an interacting quintessence model which 
shares both the same expansion history and growth history with the modified gravity model. In Sec. IIII1 
following this prescription, as an example, we explicitly demonstrate the interacting quintessence model 
which is indistinguishable with the DGP model in this sense. A brief discussion is given in Sec. IIVI in 
which some suggestions to break this degeneracy are also discussed. 



II. GENERAL FORMALISM 

Quintessence [30L l3ll is a well-known dark energy candidate. In this section, we consider the interacting 
quintessence model |32i l33l . l3~il HE HI, H3] and show how it can share both the same expansion history 
and growth history with modified gravity models, without invoking non-trivial dark energy clustering. 
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We consider a flat Friedmann-Robertson- Walker (FRW) universe. As is well known, the quintessence 
is described by a canonical scalar field with a Lagrangian density C$ = \ (d^cj)) 2 — V{<j>). Assuming the 
scalar field <f> is homogeneous, one obtains the pressure and energy density for quintessence 

P<t> = \j> 2 -v{<j>), p+ = \p + vy>), (i) 

where a dot denotes the derivative with respect to cosmic time t. The Friedmann equation reads 

3H 2 = n 2 (p m + P4> ) , (2) 

where H = a/a is the Hubble parameter; a = (1 + z)^ 1 is the scale factor (we have set ao — 1; the 
subscript "0" indicates the present value of corresponding quantity; z is the redshift); p m is the energy 
density of cold dark matter (we assume the baryon component to be negligible); k 2 = 8ttG. We assume 
that quintessence and cold dark matter interact through (32l . HH, l3~ij 

p m + 3Hp m = -KQp m (j), (3) 
P4> + 3H ( P4> + p<f,) = nQp m (j), (4) 

which preserves the total energy conservation equation p to t + 3H (p to t + Ptot) = 0. The dimensionless 
coupling coefficient Q — Q{4>) is an arbitrary function of </>. Eq. ^ is equivalent to 

<f> + 3H(j) + — = nQp m . (5) 
dip 

Using Eqs. @ — (H]), one can obtain the Raychaudhuri equation 



H = {p m + P4, + p$) = [pm + <t> J ■ (6) 

It is worth noting that due to the non- vanishing interaction, p m does not scale as a~ 3 . The above equations 
are associated with the expansion history. On the side of growth history, as shown in [33[ , the perturbation 
equation in the sub-horizon regime is 

where i5 = Sp rn /p rn is the linear matter density contrast; fl m = K 2 p m /(3H 2 ) is the fractional energy 
density of cold dark matter; and a prime denotes a derivative with respect to In a. Note that in [33T | the 
absence of anisotropic stress has been assumed, namely, in longitudinal (conformal Newtonian) gauge the 
metric perturbations $ = \&. Obviously, when Q = 0, Eq. ([7]) reduces to the standard form in general 
relativity 0, M, M M M M H 

S + 2HS = 4nGp m S. (8) 

In fact, Eq. from [33[ is valid for any Q — Q{4>) and generalizes the one of [34| which is only valid 
for constant Q. On the other hand, in modified gravity, the perturbation equation ([8|) has been modified 

to [SlilSlili^ 

I + 2H5 = ^G cSPm ~5, (9) 

where the quantities in modified gravity are labeled by a tilde "~" ; G e s is the effective local gravitational 
"constant" measured by Cavendish- type experiment, which is time-dependent. In general, G e ff can be 
written as 

Geff = G (l + j^j , (10) 
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where (3 is determined once we specify the modified gravity theory. Eq. @ can be rewritten as 

' ^ + Ts) & m 5. (11) 




2 V 3/3 

Now, we require that the interacting quintessence model shares both the same expansion history and 
growth history with the modified gravity model. That is, we identify 

H = H and 6 = 5. (12) 

Comparing Eq. ((7]) with Eq. (fTTjl . we find that 



b's' = Is 



3/3 

Note that f2 m ^ fl m in general. From Eq. ((6]), we have 



(13) 



( K( f>'f = -30 m - 2^-. (14) 
rl 



H' 



We can recast Eq. (J3j) to 



It turns out 



<>'„ = - |3 + 2^ + K Q0' j <>,„. fir,. 



K Q0' = -3 - 2^ - ^. (16) 



From Eqs. ([Till and (fT6|) . we obtain 



o, H ' -J- 

W) 2 y i+ H + n 
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Noting that <5 = 5, we can find 5 in Eq. (fi"3"|) from Eq. (jTTJ) . Once f2 m , /3, i/, and corresponding S in 
the modified gravity are given, substituting Eqs. (fTBj) and (fT7|) into Eq. (fTBf and noting Eq. (fT2|) . we 
obtain a differential equation for fi m with respect to In a. After we find f2 m (lna) from this differential 
equation, ncp'(hia) can be obtained from Eq. (fT4"|) , while H = H. Then, by using Eqs. (|T^)l and (JT^J) , 
Q(lna) = (kQ^') / (K<p') is in hand. From Eq. ©, fi = ^p^/^H 2 ) = 1 - O m . Noting Eq. JJJ and using 
Eq. (fT4|) . we find the dimensionless potential of quintessence 



fljT~ 3B I ~ + 3 



U=^=3E>[l-^ + ~^), US. 



where £7 = H/Hq. Notice that H' /H = E'/E. By integrating K0'(lna), we find K(f> as a function of In a. 
Therefore, we can finally obtain Q, U as functions of K(f>. 

In fact, this general approach proposed here is just a normal reconstruction method. For any given 
modified gravity model, we can always construct an interacting quintessence model in the framework of 
general relativity which shares both the same expansion history and growth history with this given modi- 
fied gravity model. As is well known, an interacting quintessence model can be completely described by its 
potential V((f) and the coupling Q((j>). For a given modified gravity model, its Q m , P, and H are known, 
whereas its corresponding S can be obtained from Eq. (jTTJ) . Following the procedure described in the text 
below Eq. (|17|) . one can easily reconstruct the dimensionless potential U(k<P) [which is equivalent to the 
potential V(<f>) obviously] and the coupling Q{n<p) [which is Q{<f)) in fact] for the interacting quintessence 
model. Up until now, the desired interacting quintessence model which shares both the same expan- 
sion history and growth history with the given modified gravity model has been constructed. Therefore, 
the cosmological observations might be unable to distinguish between them, unless other complementary 
probes beyond the ones of cosmic expansion history and growth history are used. 




FIG. 1: 8 — S, Q m , K(f>' , Q, U = k 2 V/Hq, and K(f> as functions of In a. See text for details. 
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III. EXPLICIT EXAMPLE 



In this section, we give an explicit example following the prescription proposed in Sec. [TTJ Here, we 
consider the DGP braneworld model (see also e.g. [H, HH, [H, S3), which is the simplest modified 
gravity model. Assuming the flatness of our universe, in the DGP model (here we only consider the 
self-accelerating branch), E = H/Hq is given by fl8l I2H |22| 



e = \jn m0 (i + zf + ci rc + yn ra = yn m0 e- 3ina + ct rc + yn ra , (19) 

where Q ra is constant. E(z = 0) = 1 requires 

ft m0 = 1 - 2y^~. (20) 

Therefore, the DGP model has only one independent model parameter. The fractional energy density of 
matter in the DGP model reads 

- _ n m0 (i + z ) 3 = n m0 e- 3lna 

E 2 (z) E 2 {\na) ' 

In addition, the (5 in Eq. (JTUJ) for the flat DGP model is given by [H, SJ H3| 

i + n 2 

/9= _i±4>. (22) 

1 - o 2 

Fitting the DGP model to the 192 SNIa data compiled by Davis el al. [lfj which are joint data from 
ESSENCE 0| and Gold07 0, we find that the best-fit parameter (l rc = 0.170 with xLin = 196.128 while 
the corresponding f2 m0 is given by Eq. (|2T)|) . Substituting this f2 r „o into Eqs. (|2"Tj) , (f2"2"|) and (fTT)|) , fi m , /3, 
and E as functions of In a are known (since our main aim is to demonstrate the prescription proposed 
in Sec. [JTJ we do not consider the errors, even one can use any fi TO o here for demonstration). Noting 
that H' /H = E'/E, following the prescription proposed in Sec. HH we can easily construct the desired 
interacting quintessence model which shares both the same expansion history and growth history with the 
corresponding DGP model. At the first step, we obtain 5 = 6 from Eq. (fTTj) . As is well known, 8 = 6 = a at 
z > 1 (see e.g. 0IH). Thus : 

we use the initial condition 5' = 5 = ai n i at Zi n i = 1000 for the differential 
equation (llip . The resulting 6 = S as a function of In a is shown in Fig.[T] At the second step, substituting 
Eqs. UH) and (TT7| into Eq. (TTS)) . we obtain Q, m as a function of In a from the resulting differential equation. 
Note that cold dark matter can decay to quintessence through interaction, f2 m is unnecessary to be 1 at 
high redshift. So, for demonstration, we choose the initial condition f2 m (z = Zini) = 0.995 at z; m i = 1000 
for the differential equation of fl rn . Different values of £l m (z — Zi n i) only mean different displacements of 
Kip' and the dimensionless potential U = k 2 V/Hq at Zj n j. The resulting £l m as a function of In a is also 
shown in Fig. [TJ Then, following the prescription proposed in Sec. HH it is straightforward to obtain K<p' , 
Q, U = k 2 V/Hq, and n<j> as functions of In a, while for demonstration we choose the negative branch for 
K<fi' , and choose 4>o = when we get K(p. They are also shown in Fig. [TJ Once we obtain Q, U and n<j) as 
functions of In a, it is easy to find Q and U as functions of K<f>. The results are shown in Fig. [5] 

In short, here we have faithfully followed the reconstruction method proposed in Sec.[TTj and successfully 
constructed an interacting quintessence model which shares both the same expansion history and growth 
history with the DGP model whose single model parameter (l rc = 0.170 which is the best- fit to the 
192 SNIa data compiled by Davis et al. [10| for example. The potential V{(j>) and coupling Q(4>) which 
are required to describe the reconstructed interacting quintessence model have been presented in Fig. [5] 
through the equivalent U(n(f>) and Q(n<p). There are no analytical expressions for them and instead the 
reconstructed U(n(f>), Q(n(f>) are given in numerical forms. 



IV. CONCLUSION AND DISCUSSION 



In summary, we proposed a general approach to build an interacting quintessence model which shares 
both the same expansion history and growth history with the modified gravity model. Therefore, the 
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cosmological observations might be unable to distinguish between them, unless other complementary 
probes beyond the ones of cosmic expansion history and growth history are used. As an example, we 
also explicitly demonstrated the interacting quintessence model which is indistinguishable with the DGP 
model in this sense. In fact, this proposed prescription is also valid for other modified gravity models, 
such as f(R) theories, braneworld-type models, scalar-tensor theories (including Brans-Dicke theory), 
TeVeS/MOND models, and so on [43|, HE] • Of course, one can also extend the interacting quintessence 
model to other interacting dark energy models. 




0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 

K(p K<f> 



FIG. 2: Q and U = k 2 V/Hq as functions of K(j>. 



In this note, without invoking non-trivial dark energy clustering (e.g. non- vanishing anisotropic stress), 
we find that the interaction between dark energy and dark matter could also make dark energy model 
and modified gravity model indistinguishable. How can we break this degeneracy? Firstly, we should 
carefully check the evidences of the interaction between dark energy and dark matter [H, \4^. If 
this interaction does not exist, the combined probes of cosmic expansion history and growth history 
might be enough to distinguish between dark energy and modified gravity. It is worth noting that the 
current observational constraints on the interaction between dark energy and dark matter in the literature 
(e.g. [H, [2i| El, 46]) considered other interaction forms (e.g. oc Hp m , Hp c i e or Hp to t) which are different 
from the one of the present work [cf. Eqs. ^ and Q], and hence cannot be used to compare with the 
interaction considered here. Therefore, it is of interest to consider the observational constraints on the type 
of interactions oc Q{4>)p m 4> in future works. Secondly, to break the degeneracy, the other complementary 
probes beyond the ones of cosmic expansion history and growth history are desirable. For instance, these 
complementary probes might include local tests of gravity, high energy phenomenology, and non-linear 
structure formation. Thirdly, in addition to the linear matter density contrast S(z), one can also test 
the metric perturbations <I> and 'F by using the relationship between gravitational lensing and matter 
overdensity [5(J. For the interacting quintessence model $ = *F, whereas for the DGP model $ ^ >F. 
Thus, it is possible to distinguish between them by using delicate measurements. We consider that this 
issue deserves further investigations and believe that a promising future is awaiting us. 
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